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Abstract 


There  is  a  requirement  to  soft  catch,  wititout  exceeding  a  deceleration  rate  of  1,600  g, 
and  within  less  than  300  m,  an  unmodified  102-lb  SAD  ARM  (Sense  and  Destroy  Armor) 
projectile  fired  at  840  m/s.  This  report  presents  a  soft-recovery  concept  and  its 
numerical  simulation.  The  concept  entails  aerodynamic  deceleration  of  the  projectile  in 
a  long  tube  attached  to  the  gim  barrel.  The  midsection  of  the  tube  is  boimd  between  a 
diaphragm  and  a  free  piston  and  is  prepressurized  to  about  2  MPa.  As  the  projectile 
enters  the  tube,  die  shock  wave  preceding  it  rupttires  the  diaphragm  and  the  projectile 
decelerates  as  high  pressxue  builds  between  it  and  the  free  piston.  The  piston 
disengages  and  travels  forward  scooping  water.  The  waterlog  that  forms  in  front  of  the 
piston  effectively  increases  the  piston's  mass  and  also  induces  braking  force  because  of 
the  water  friction  with  the  tube  wall.  The  projectile's  deceleration  is  controlled,  and 
eventually  the  projectile  exits  the  tube  with  a  velocity  of  10  m/s.  The  numerical 
simulation,  based  on  the  method  of  characteristics,  incorporates  unsteady 
one-dimensional  fluid  dynamics  titiat  captures  the  extensive  wave  dynamics.  This 
report  details  the  effects  on  the  projectile's  deceleration  of  die  midsection  length,  initial 
pressure,  and  the  water  mass.  From  the  simulation,  it  is  possible  to  soft  capture  the 
SAD  ARM  projectile  within  120  m. 
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1.  Introduction 


1.1  The  Requirement  for  a  Soft  Recovery  System 

With  advances  in  microelectronic  and  electro/optic  sensors,  smart  projectiles 
have  become  a  reality.  Such  a  projectile  is  the  U.S.  Army's  XM898  SADARM 
(Sense  and  Destroy  Armor)  155  mm.  Because  of  the  sensitive  components  and 
packaging  inherent  to  these  designs,  projectile  tests  to  evaluate  performance 
and/or  failure  analysis  are  very  expensive.  Owing  to  the  projectile's  complexity 
and  thin-wall  construction,  a  dummy  or  malfunctioned  projectile  recovered  after 
being  fired  downrange  may  be  too  damaged  for  useful  analysis.  Hence  arose  the 
need  for  a  soft  recovery  system  (SRS)  for  the  XM898  and  similar  projectiles.  Soft 
recovery  is  defined  here  as  the  recovery  of  the  projectile  in  a  manner  that  the 
projectile  does  not  exceed  certain  deceleration  limits  and  damage  thresholds. 
The  foremost  proponents  of  the  SRS  have  been  Ami  Frydman  [1]  and  Donald 
Carlucci  [2].  The  major  specifications  for  the  SRS  are  (1)  the  projectile  is  to  be 
fired  from  a  standard  gun  barrel  with  standard  propellant  charge,  (2)  the 
projectile  shall  not  be  modified  in  any  way,  (3)  the  SRS  shall  be  able  to  soffly 
catch  a  projectile  fired  at  840  m/s  within  200  m  from  the  muzzle  without 
exceeding  a  deceleration  level  of  1,600  g,  and  (4)  the  SRS  shall  have  a 
rapid-tumaroxmd.  The  utilization  of  the  SRS  will  enable  tiie  verification  of 
latmch  integrity  and  functionality  of  the  projectile's  components,  by  measuring 
the  functional  eind/  or  structural  performance  of  the  components  during  the  first 
1,200  cal.  of  the  projectile's  travel.  The  principal  author  of  this  report  has 
proposed  a  novel  SRS  concept  that  fulfills  the  SRS  specifications  and  is  believed 
to  have  important  advantages  when  compared  to  other  concepts.  This  report 
details  a  numerical  model  and  its  predictions  for  the  flow  dynamics  in  the 
proposed  SRS — work  funded  by  the  SADARM  Program  Management  Office  to 
verify  the  feasibility  of  the  concept  and  provide  a  tool  for  its  engineering. 

1.2  Soft  Recovery  Concepts  and  Systems 

Basically,  all  of  the  techniques  for  the  soft  recovery  of  fired  projectiles  had  been 
established  over  20  years  ago.  Some  of  the  techniques  were  actually  tried  out; 
others  existed  only  in  the  form  of  patents.  However,  only  a  few  of  the  techniques 
comply  with  the  present  specifications  for  the  SRS.  Wright  [3]  gave  a 
comprehensive  list  of  techniques  that  existed  up  to  the  early  1970s.  Paul  Baer  [4] 
of  the  U.S.  Army  Ballistic  Research  Laboratory  (BRL),*  Aberdeeii  Proving 
Ground,  MD  proposed  an  SRS  for  the  155-mm  carmon.  He  suggested  firing  the 


*  The  U.S.  Army  Ballistic  Research  Laboratory  (BRL)  was  deactivated  30  September  1992  arvd 
subsequently  became  part  of  the  U.S.  Army  Research  Laboratory  (ARL)  on  1  October  199Z 
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projectile  in  a  side-vented  gun  such  that  the  projectile  exits  with  a  low  velocity 
into  a  recovery  tube  filled  with  compressed  air.  The  rifled  gun  tube  was  to  be 
connected  to  the  recovery  tube  via  a  "rotating  band  squeezer"  followed  by  a 
diaphragm  that  would  retain  the  compressed  air.  Paul  Baer  simulated 
numerically  the  flow  in  his  proposed  SRS,  but  used  a  lumped  parameter 
approach  that  completely  ignored  the  very  important  wave  motion  that  can 
produce  large  peaks  of  projectile  deceleration.  A  major  drawback  of  Baer's  SRS 
is  the  requirement  to  modify  the  gun  by  the  side  vent.  In  the  early  1980s, 
Honeywell  Inc.  [5]  constructed  an  SRS  for  BRL  and  the  system  was  tested 
successfully  with  a  155-mm  cannon  [6].  Although  the  Honeywell  system  had  the 
advantage  of  being  able  to  accommodate  multiple  calibers,  it  had  a  major 
drawback -it  required  a  modification  of  the  projectile.  The  system  employed  a 
water  scoop  mounted  to  the  projectile  that  slowed  the  projectile  via  momentum 
exchange  with  water  in  a  trough.  The  projectile  also  was  modified  during  flight 
because  the  technique  required  the  stripping  of  the  projectile's  A  (rotating)  band. 
A  physically  sotmd  concept  for  an  SRS  is  the  Ballistic  Compression  Decelerator 
patented  by  the  McDoimell  Douglas  Corporation  in  1972  [7].  The  projectile  is 
fired  into  a  prepressurized  succession  of  "decelerator  tubes"  that  are  separated 
by  multiple  diaphragms  such  that  the  pressure  builds  up  ahead  of  the  projectile 
(tiius  decelerating  it),  and  the  diaphragms  rupture  before  being  pierced  by  the 
projectile. 

An  SRS  based  on  the  McDonnell  Douglas  patent  is  presently  in  operation  [1, 2]  at 
the  German  firm  RheinmetaU  W&M,  and  this  SRS  is  ihe  main  competitor  for  the 
SRS  proposed  in  this  report.  It  can  stop  an  840  m/s  155-mm  projectile  in  200  m. 
The  main  disadvantage  of  the  mtdtiple  decelerator  tube  concept  is  that  it 
employs  long  tubes,  and  mxiltiple  diaphragms.  Therefore,  it  is  expensive  to 
construct  and  it  occupies  a  large  amount  of  real  estate.  A  variation  on  the 
ballistic  compression  concept  is  the  McDonnell  Douglas  1976  patent  [8]  titled 
"Projectile  Recovery  System  With  Quick  Opening  Valves."  Here,  quick  opening 
valves  that  unseal  the  tube  as  the  projectile  approaches  them  replace  the 
diaphragms.  This  concept  is  considered  impractical  because  of  its  reliability. 
Two  patents  that  do  not  rely  on  ballistic  compression  are  the  "Rifled  Soft 
Recovery  System,"  [9]  and  "Multicaliber  Projectile  Soft  Recovery  System"  [10]. 
The  rifled  concept  relies  on  a  spinning  projectile  being  accurately  reengaged  in  a 
recoiled  rifled  tube.  This  concept  requires  great  precision  of  the  firing  cycle  and 
tube  alignment;  it  is  not  proven  and  will  probably  subject  the  projectile  to 
balloting  and  excessive  stresses.  The  multicaliber  concept  has  the  advantage 
(unlike  ballistic  compression)  of  handling  multi  calibers  for  the  same  SRS.  A 
deformable  element  is  placed  in  the  path  of  the  projectile  and  the  projectile 
becomes  embedded  in  the  element.  Mechamcal  forces  are  applied  to  stop  the 
element.  This  concept  lacks  the  deceleration  control  that  the  ballistic 
compression  technique  has,  and  it  subjects  the  projectile  to  excessive  forces  upon 
embedding. 
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2.  The  Proposed  Soft  Recovery  System 


The  proposed  SRS  concept  for  the  155-iiim  projectile  is  based  on  the  ballistic 
compression  principle.  It  is  sketched  in  Figure  1.  Unique  to  the  concept  is  the 
incorporation  of  a  water-controlled  free  piston  in  the  decelerator  tube.  The 
principle  of  operation  is  as  follows.  The  atmospheric  transition  tube  (ATT)  has 
rifling  that  continues  the  barrel's  rifling.  The  rifling  is  tapered  and  it  transitions 
to  a  smooth  barrel.  The  spin-stabilized  155-mm  projectile  enters  the  ATT  where 
the  tapered  rifling  in  the  tube  gradually  compresses  the  projectile's  rotating  band 
such  that  the  band  engraving  is  squeezed  out  and  a  good  seal  is  formed  between 
the  band  and  the  smooth  wall  portion  of  the  ATT.  The  tapered  rifling  is  believed 
to  alleviate  balloting  of  the  projectile  in  the  tube.  The  band  sealing  assures  that 
compressed  gas  in  front  of  the  projectile  will  not  leak  around  the  projectile.  The 
combustion  gases  behind  the  projectile  escape  through  the  vent  holes  in  the  ATT 
and  the  pressure  on  the  projectile's  rear  face  effectively  drops  to  atmospheric 
value.  The  shock  wave  that  precedes  the  projectile  ruptures  the  750-psi 
diaphragm  that  holds  the  pressurized  air  and  enters  into  the  pressurized  tube. 
When  the  shock  wave  reaches  the  free  piston  and  reflects  from  it,  the  free  piston 
overcomes  the  shear  ring  that  holds  it  against  the  pressurized  air.  The  free  piston 
then  moves  down  the  guide  tube  pushing  ahead  a  growing  column  of  water. 
The  free  piston  is  made  of  tight  material,  e.g.,  plastic  such  as  high-density 
polyethylene.  The  free  piston  has  to  be  tight  enough  so  that  initially  it  accelerates 
quickly  and  thus  prevents  high-pressure  spikes  from  developing  in  the 
pressurized  tube  section.  The  water  effectively  increases  the  free  piston  mass 
with  time  and  also  adds  frictional  force  due  to  the  water  boundary  layer  on  the 
tube  wall.  The  growth  of  the  free  piston  mass  and  the  added  water  friction 
retard  the  free  piston  motion  and  thus  regulate  the  air  pressure  between  the  free 
piston  and  the  projectile  to  a  sustained  level  between  2,000  and  4,000  psi.  This 
pressure  is  enough  to  slow  the  projectile  from  840  m/s  to  10  m/ s  within  100  m 
without  exceeding  the  1600-g  deceleration  limit.  The  water  and  compressed  air 
escape  in  the  slotted  brake  tube,  and  the  projectile  is  stopped  in  the  tapered 
lining  of  the  brake  tube.  In  operation  of  the  SRS,  the  initial  pressure  and  water 
mass  are  easily  adjusted  to  allow  fine  control  of  the  projectile  exit  velocity  into 
the  brake  tube. 

The  expendable  parts  in  this  SRS  concept  are  the  diaphragm,  the  shear  ring,  and 
the  free  piston.  Because  the  proposed  SRS  length  is  much  shorter  than  the 
aforementioned  RheinmetalTs  SRS  and  it  employs  fewer  expendable  parts,  it  is 
believed  that  the  construction  and  operation  of  the  proposed  SRS  is  significantly 
cheaper  than  Rheinmetall's  SRS. 
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3.  Mathematical  Model  of  the  Flow  Dynamics  in  the  SRS 
Tube  Sections 


A  numerical  model  was  constructed  to  simulate  the  transient  flow  fields 
generated  in  the  SRS  tube.  The  model  consists  of  four  separate  zones  linked  by 
botmdary  conditions  across  the  projectile,  the  free  piston,  the  two  shock  waves, 
and  the  exit  plane  (Figure  2). 


traveling 

projectile 


□ 


Free  Piston  & 
Water  Column 

-► 


traveling  shock  wave 


distance  between  Traveling  Shock  Wave 
A  and  B  grows  with 
time 


Figure  2.  The  zonal  construction  of  the  SRS  tube  flow. 


Zone  1  is  between  the  projectile  and  the  traveling  shock  wave  that  originally  is 
generated  by  the  supersonic  entry  of  the  projectile  into  the  tube.  Zone  2  is  the 
region  between  this  traveling  wave  and  the  free  piston.  Initially,  the  shock  wave 
leading  the  projectile  reflects  from  the  rupture  diaphragm  (the  diaphragm  is 
located  in  zone  2,  although  it  is  not  shown  in  Figure  2).  The  diaphragm  should 
rupture  upon  reflection,  thus  initiating  the  transient  flow  in  the  pressurized  gas 
between  the  diaphragm  and  the  free  piston.  When  the  pressure  on  the  free 
piston  exceeds  a  certain  value,  the  free  piston  overcomes  the  shear  ring  that 
holds  it  against  the  pressurized  air  and  begins  its  motion.  The  traveling  wave 
may  reflect  a  few  times  between  the  free  piston  and  the  projectile  imtil  the  free 
piston  has  exited  the  tube.  In  the  numerical  solution  procedure,  the  decision  to 
reflect  the  wave  from  either  the  projectile  or  the  free  piston  is  made  when  the 
thickness  of  zones  1  or  2  becomes  less  than  1%  of  the  distance  between  the 
projectile  and  free  piston. 

Zone  3  is  between  the  front  of  the  water  column  and  the  traveling  shock  wave 
generated  by  the  fast  acceleration  of  the  free  piston.  The  water  column  thickens 
as  the  free  piston  sweeps  the  water  layer  ahead  of  it.  The  water  column  forms  in 
front  of  the  piston  because  the  piston  velocity  is  1  order  of  magnitude  greater 
than  the  water  gravity  wave  (for  wave  amplitude  that  equals  the  tube  diameter). 
Zone  4  is  between  the  traveling  wave  and  the  tube's  exit.  Zone  4  disappears 
when  the  traveling  wave  reaches  the  exit,  and  at  this  time  zone  3  becomes  the 
region  between  the  free  piston  and  the  exit.  Then,  zone  3  disappears  when  the 
front  of  the  water  column  reaches  the  exit  plane  and  the  free  piston  exits  the 
tube.  Zone  2  disappears  when  the  traveling  wave  between  zones  2  and  1  reaches 
the  exit,  and  at  this  time  zone  1  becomes  the  region  between  the  projectile  and 
the  exit. 
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The  time-dependent  governing  equations  in  each  zone  assume  isentropic  flow  of 
a  Noble- Abel  gas  with  constant  molecular  weight,  covolume,  and  ratio  of  specific 
heats.  Furthermore,  the  flow  is  assumed  to  be  one-dimensional  and  inviscid, 
which  of  course  neglects  heat  transfer  to  the  tube  wall.  The  one-dimensional 
treatment  of  the  problem  necessarily  assigns  a  flat  front  to  the  projectile,  a  front 
that  in  reaUty  is  an  elongated  ogive.  Obviously,  boundary  layer  effects  are  not 
accoimted  for.  In  practice,  the  growth  of  boundary  layers  in  the  long  sections  of 
the  tube  will  constrict  the  flow  and  modify  tine  strength  of  the  traveling  shock 
waves.  Intuitively,  viscosity  and  boundary  layer  effects  will  tend  to  dissipate  the 
wave  motion  in  the  tube  thus  resulting  in  lower  peak  pressure  loading  (and 
hence  deceleration)  on  the  projectile.  Thus,  the  simplified  assumptions  are  likely 
to  result  in  conservative  peak  pressure  values.  The  present  model  incorporates 
the  most  important  aspects  of  the  flow  physics  and  it  is  therefore  adequate  for 
ihe  purpose  of  feasibility  studies  of  the  SRS  concept,  and  for  a  parametric  study 
of  the  concept. 

With  pressure  P,  density  p,  velocity  u  (with  respect  to  laboratory  coordinates), 
and  sound  speed  a,  these  equations  can  be  written  in  characteristic  form  as 


dP  A 

—  ±p-a - =  0, 

dt  dt 


(1) 


which  are  integrated,  respectively,  along  the  right  (+)  and  left  (-)  nmning 
directions. 


dZ  TTJ. 

—  =  U±a, 
dt 


(2) 


where  U  is  the  material  velocity  relative  to  the  coordinate  velocity  Vc,  i.e.,  U  -  u  - 
Vc.  The  energy  equation  can  be  written  in  the  form 


dt' 


(3) 


which  is  integrated  along  the  streamline 


(4) 


With  the  temperature  T,  covolume  b,  ratio  of  specific  heats  y  (=Cp/  Cv),  and  gas 
constant  SR  (=  tmiversal  gas  constant/molecular  weight),  the  Noble-Abel 
equation  of  state  is  simply 


P  SRT 
p  1-p-b' 
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which  dictates  the  sound  speed  a, 


2  dP 
a  = 


dp 


entropy 


y-P/p 

1-pb 


(6) 


There  are  three  independent  variables  here:  u,  P,  and  p.  Two  equations  are 
provided  by  the  compatibility  conditions  in  equation  1  along  the  right  and  left 
characteristic  directions  given  by  equation  2.  The  third  is  the  energy  equation 
given  in  equation  3  which  can  be  solved  along  the  streamline  in  equation  4.  The 
temperature  T  follows  directly  from  equation  5. 

As  previously  mentioned,  when  the  coordinate  system  is  in  motion,  the 
characteristic  directions  must  be  calculated  with  the  velocity  relative  to  the 
moving  coordinate,  i.e.,  U  =  u  -  Vc.  For  all  calculations  discussed  here,  the  end 
points  of  each  coordinate  zone  remain  attached  to  the  adjacent  moving 
boimdary,  and  the  local  grid  points  within  the  zone  are  distributed  tmiformly 
between  the  instantaneous  boundary  locations  (like  an  "accordion").  Hence,  the 
local  coordinate  velocity  varies  linearly  between  the  two  boundary  velocities. 
This  is  illustrated  schematically  in  Figure  3  for  zone  1  which  is  attached  to  the 
projectile  front  face  and  the  traveling  shock  wave. 


U,(l+dt) 


time  s  t+dt 


u.w(t+dt) 


1  2  3  4  5  6 


□ 


“pW . 


time  =  t 


1  2  3  4  5  6 


Projectile  Traveling  Shock  Wave 


Figure  3.  Illustration  of  numerical  coordinate  motion  with  time. 

The  numerical  solution  of  equations  1  and  3  is  determined  with  the  method-of- 
characteristics  (MOC)  which  to  a  great  extent  will  faithfully  reproduce  local 
wave  motion  without  spurious  numerical  oscillations.  The  MOC  routines 
employed  here  were  taken  from  Kooker  [11]  where  all  interpolation  is  done  with 
Akima  splines  [12].  All  boimdary  values  are  also  obtained  with  special  routines 
based  on  MOC.  Motion  of  the  projectile  and  the  free  piston  within  the  tube  is 
assumed  to  be  frictionless. 
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3.1  Solution  Methodology 

In  general,  the  flow  is  solved  sequentially  from  zones  1  to  4.  For  each  zone,  the 
boundary  conditions  must  be  solved  first  before  the  interior  flow  solution.  The 
solution  is  known  at  time  t=0  and  it  is  marched  from  time  t  to  time  t+dt  using  the 
characteristic/free-stream  construct  shown  in  Figures  4-8.  All  notations  are 
noted  in  these  figures. 

3.1.1  The  Projectile  Boundary  Conditions  (Figure  4) 


Figure  4.  Construct  of  the  projectile's  boundary. 


From  equations  1  and  2,  upstream  along  the  4-1  left-running  characteristic: 

P,-P4-(p-a)4i -(01-04)  =  0, 


(7) 


where 


(p-a)4+(p-a)i] 


From  Newton' s  Second  Law: 


1*20  + 1*2 


Yr  d(M-v) 

^  dt 


l-A„  =M„  • 

0 

1 

IP  P  1 

1  dt  J 

(8) 


8 


71  * 

Mp  and  Ap  =  — - —  are  the  projectile  mass  and  cross-section  area. 

From  equations  3  and  4,  along  the  streamline  10-1; 

Pi-Pio-am(pi-Pio)  =  0/  (9) 

fa^+aM 

where  afj^  =  <  — - —  >  and  a  is  given  by  equation  6. 

2i 

> 


The  pressures  behind  the  projectile  P20  and  P2  are  given  (assumed  to  be 
atmospheric).  The  values  at  point  10  are  known  from  the  previous  step,  and  the 
values  at  point  4  are  found  using  Akima  splines  interpolation  for  this  point. 
Thus,  the  unknowns  in  equations  7-10  are  Uj ,  Pj ,  and  pj ,  and  they  are  fotmd  by 
guessing  an  initial  Pj  value  and  solving  equations  7-9  iteratively. 


3.1.2  The  Free-Piston  Boundary  Conditions  (Figure  5) 


Figure  5.  Construct  of  the  free-piston  boundaries. 

The  solution  methodology  is  similar  to  that  of  the  projectile,  with  the  exception 
that  the  right  running  characteristic  has  to  be  solved,  and  the  water  column  mass 
and  friction  have  to  be  accounted  for  in  the  force  balance  equation. 
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nA 


From  equations  1  and  2,  downstream  along  tiie  2-3  right-running  characteristic: 

P2-P3+(p.a)32-(u2-U3)  =  0,  (10) 

where 

(p*a)32  = 

From  equations  3  and  4,  along  the  streamline  20-2: 

1*2  ~^20  "^220  ■(P2  “P2o)“ 

where 


j'(p-a)3  +  (p-a)2l 


Similarly,  along  the  streamline  10-1: 

Fi  ~Pio  ~^iio  '(Pi  “Pio)" 

where 


Upstream  along  tiie  7-1  left-nmning  characteristics: 

P,-P2-(p-a)7i-(u,-U7)  =  0, 

where 

From  Newton' s  Second  Law: 


(12) 


(13) 


■  (Lw  -u’  fl..}- A„  +M*,)-u, 

—  (Mpp  +  M^jo)  •  Uio}/ 

where  M^io  and  M^,  are  the  masses  of  the  water  column  at  time  t  and  t+dt, 
is  the  water  column  length,  and  Mpp  and  App  =  Ap  are  the  free  piston  mass  and 
cross-  section  area. 
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The  frictional  force  due  to  the  water  column. 


where  the  friction  coefficient  is  taken  from  [13]  for  turbulent  flow  in  a  tube  (the 
Prandtl  universal  law  for  smootiti  circular  tube). 

f  = - - r - ^  where  Re  =  ^  is  the  Reynolds  number. 

16-(log,o(ReVf)-0.4j  Pw 

There  are  six  unknowns  in  equations  10-14,  namely:  Uj ,  Pj ,  pj ,  and  U2 ,  P2 ,  P2  / 
where,  of  course,  sotmd  speeds  ai  and  a.z  are  found  from  equation  6.  The 
remaining  equation  relates  u,to  U2,  i.e.,  Ui=  (1  +  K)  •U2,  where  the  factor  K  is 
the  initial  ratio  of  water  to  air  volumes  in  the  guide  tube.  This  solution  requires 
iteration  similar  to  the  projectile  solution. 

3.1.3  Solution  for  a  Grid  Point  in  the  Inner  Flow  (Figure  6) 


A  POINT  IN  THE  INNER  FLOW 


dt 


/I 

*  • 

'  I 


U  +a 


/ 

/U 


U  -a 


Characteristic 

Streamline 


z 

-> 


Figure  6.  Construct  of  an  inner  flow  point. 


The  timer  flow  solution  is  based  on  two  characteristics  lines  (right-running  and 
left-rtmning)  and  one  streamline.  The  unknowns  are  Uj ,  p, ,  and  Pj  which  are 
found  from  equations  similar  to  equations  10,  11,  emd  13  (where,  again,  the 
sound  speed  is  found  from  equation  6). 
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3.1.4  The  Shock  Wave  Boundary  Jump  Conditions  (Figure  7) 


Figure  7.  Construct  of  the  traveling  shock  waves. 

Because  of  tihe  projectile's  supersonic  velocity,  a  shock  wave  leads  the  projectile's 
motion.  Once  the  projectile  enters  the  catch  tube  (that  has  the  same  diameter  as 
the  projectile),  the  shock  wave  detaches  from  the  projectile  and  becomes  a 
traveling  wave.  A  traveling  shock  wave  will  also  originate  in  front  of  the  free 
piston  as  the  piston  accelerates  to  supersonic  velocity. 

The  values  at  point  1  (Figure  7)  are  uniquely  determined  by  the  upstream  flow 
field  and  are  fotmd  identically  to  the  inner-point  solution.  The  remaining 
tmknowns  are  Uj,  P2,  P2/  and  (the  shock  wave  velocity),  requiring  four 
equations  for  a  solution.  The  characteristic  line  trailing  the  shock  wave  on  the 
high-pressure  side  provides  one  equation.  The  remaining  three  equations  are 
jump  conditions  across  die  discontinuity  to  conserve  mass,  momentum,  and 
energy. 
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P2  +  p2U2-(u2-Usw)  =  Pl+PlUr(Ui-Usw);  (16) 


Energy; 

e2+^  +  i-(u2-Usw)^  =e, +-^  +  ~(u2-Usw)^  (17) 

P2  2  p,  2 

where  e  =  Cy  •  (T  -  To)  and  T  is  related  to  P  and  p  through  the  equation  of  state 
(equation  5). 

Equations  15-17  can  be  rearranged  such  that  together  with  die  characteristic 
equation,  they  result  in  two  nonlinear  algebraic  equations  that  include  only 
and  P2  as  unknowns.  The  numerical  solution  is  obtained  by  iteration. 

3.1.5  The  Tube  Exit  Boundary  (Figure  8) 


Figure  8.  The  tube  exit  boundary. 


If  the  exit  flow  is  supersonic  or  sonic,  the  unknowns  areUj ,  Pj ,  P, ,  the  solution  is 
like  an  inner-point  solution.  If  the  flow  is  subsonic,  Pj  is  atmospheric,  and 
u,  and  p,  are  found  from  the  solution  of  the  right-running  3-1  characteristic  and 
the  streamline. 


13 


The  solution  methodology  is  as  follows.  An  initial  value  of  p,  is  guessed,  the 
equations  are  solved,  p,  is  updated,  and  the  point  4  is  found.  If  the  point  lies  left 
of  the  exit  plane,  the  flow  is  supersonic,  and  tiie  equations  are  re-solved  using  the 
two  characteristics  3-1  and  4-1  until  p,  converges.  If  the  point  4  Hes  right  of  the 
exit  plane,  the  flow  is  subsonic,  and  the  equations  are  re-solved  using  only  the 
3-1  characteristic  and  the  assumption  that  the  static  pressure  at  the  exit  plane  is 
atmospheric. 

3.1.6  The  Shock  Wave  Reflection  (Figure  9) 


Figure  9.  Flow  regions  in  a  reflecting  shock  wave. 

The  shock  wave  that  separates  zone  1  from  zone  2  will  reflect  from  the  rupture 
disc,  thereby  causing  a  pressure  differential  across  the  disc  that  is  sufficient  for 
its  rupture.  During  the  flow  process,  the  wave  also  reflects  from  both  the 
projectile's  face  and  the  free  piston's  rear. 

Figure  9  shows  schematically  the  shock  wave  reflection  process  from  a  botmdary 
that  may  have  a  velocity  u,.  The  flow  conditions  in  the  "front"  and  "behind" 
regions  are  known  from  the  preceding  shock  wave  boundary  solution.  The 
boundary  velocity,  u^,  which  is  also  the  gas  velocity  in  Ihe  reflected  region  at  the 
boundary,  is  known  from  the  projectile  or  free  piston  bormdary  solutions.  Thus, 
the  unknowrrs  are  p,,  P^,  and  Ugwr  (the  reflected  shock  wave  velocity).  These 
unknowns  are  easily  solved  using  the  conservation  of  mass,  momentum,  and 
energy  as  in  equations  15-17. 
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Mass: 

Pb'C'^b  '^SWr)”Pr '(^r  '^SWr)/ 

(18) 

Momentum: 

Pr  +Pr  -Ur  ’(Ur  -Uswr)==Pb  +Pb  ’Ub  ’(Ub  "UgWr); 

(19) 

Energy: 

+  ^  + 1  •  (Ur  -Uswr)'  =  +  ^  + 1  •  (Ub  -Ug^,)'  • 

Pr  2  Pb  2 

(20) 

3.2  The  Initial  Values  of  the  Flow  in  the  Tube 


At  time  t=0,  the  projectile  enters  the  tube  with  supersonic  velocity  (up  = 
840  m/s).  The  velocity,  of  the  resultant  shock  wave  is  obtained  by 

employing  an  impulsive  piston-motion  equation  [14] 


Ut 


Y  +  1 


M 


sw 


M 


(21) 


SW, 


where  ^  and  the  subscript  2  denotes  zone  2  (Figure  2). 

3-2 

The  jump  conditions  across  the  shock  wave  (equations  15-17)  yield  the  initial 
flow  parameters  in  zone  1. 


As  stated  before,  the  rupture  disc  is  designed  to  rupture  once  the  shock  wave 
reflects  from  it.  At  this  time,  the  flow  conditions  at  the  location  occupied  by  the 
rupture  disc  are  found  using  shock-tube  equations  [14]. 


p  _2Y2P’^nip  Y2P 
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(22) 

(23) 

(24) 
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where  the  subscript  r  refers  to  the  conditions  at  the  reflected  region  of  the 
incident  shock  wave  (Figure  9),  the  subscript  2P  refers  to  pressurized  zone 
between  the  ruptured  disc  and  the  free  piston,  and  the  subscript  rup  refers  to  the 
conditions  initiated  at  the  rupture  disc  location.  In  equations  22-24,  the 
unknowns  are  Pn,p ,  u^p ,  and  . 

“swmp  =  ^nip  ■  ^2P  the  speed  of  the  shock  wave  generated  by  rupturing  the  disc. 

In  reaUty,  at  the  time  of  the  rupturing,  the  projectile  is  very  close  to  the  rupture 
disc,  and  hence,  the  shock  wave  that  reflects  from  the  rupture  disc,  reflects 
shortly  from  the  projectile  and  overtakes/ overwhelms  the  rupture  disc  shock 
wave. 


4.  Simulation  Results 


4.1  Coding  and  Run  Parameters 

The  numerical  solution  was  coded  in  FORTRAN  77  using  the  Watcom  [15] 
Compiler.  The  code  is  run  via  an  executable  file  in  Microsoft  Windows 
(Windows  NT  4.0,  Windows  95,  or  greater)  environment.  The  program  input  is 
done  interactively,  and  output  data  is  automatically  stored  in  ASCII  files  on  the 
computer  hard  disc.  The  program  prompts  for  lengths  of  the  SRS  tube  sections; 
the  weights  of  the  projectile,  free  piston,  and  water;  the  initial  pressure  in  zone  3; 
the  projectile  velocity;  the  hold  pressures  of  the  rupture  diaphragm  and  shear 
ring.  Other  input  parameters  include  die  number  of  grid  points  per 
computational  zone,  time  constants,  rtm  options,  and  output  frequency. 
Typically,  the  number  of  grid  points  per  computational  zone  varies  from  11  to  51 
points  where  computational  zones  1  and  2  are  more  finely  resolved.  These  grid 
resolutions  were  found  to  be  adequate -higher  resolution  did  not  produce 
materially  different  results.  The  time  steps  of  the  solution  were  determined  from 
the  Courant-Friedrichs-Lewy  condition  [16].  When  the  free  piston  or  the  shock 
waves  exit  the  tube,  this  causes  a  step  pressure  drop  at  the  exit  plane,  and  it  was 
necessary  for  numerical  stability  to  represent  the  step  pressure  as  a  time-wise 
exponential  function.  Typically,  the  time  constant  for  the  exponential  function 
was  chosen  to  be  from  3  to  7  time  steps.  Run  options  include  various  default 
parameters;  and  in  addition  the  choices  of  ignoring  the  water  column  friction 
and/or  ignoring  the  transition  tube  and  the  rupture  disc.  Another  run  option  is  a 
Regula  Falsi  [17]-based  iteration  of  either  the  initial  pressure  in  zone  3,  or  the 
total  mass  of  the  distributed  water,  for  a  desired  projectile  exit  velocity.  For  the 
simulation  runs  described  here,  the  free-piston  weight  is  fixed  at  5  kg  that  is  a 
practical  weight  for  a  polyethylene-made  piston.  Also,  the  projectile  weight  and 
muzzle  velocity  are  fixed  at  the  operational  values  of  46  kg  and  840  m/  s.  The 
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pressure  rating  of  the  rupture  diaphragm  is  chosen  to  be  750  psi  (5.17  MPa)  and 
the  shear  ring  holds  up  to  500-psi  (3.45  MPa)  loading  on  the  free  piston.  The 
transition  tube  length  is  fixed  at  4.5  m. 

4.2  Typical  Flow  Dynamics  With  and  Without  Inclusion  of  the 
Rupture  Disc  Process  (Figure  10) 


(6 
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Figure  10.  Dynamics  in  the  SRS — effects  of  the  transition  tube/rupture  diaphragm 
process. 

Figure  10  demoristrates  the  flow  dynamics  in  an  SRS  where  the  transition  tube  is 
4.5  m  long  (shown  at  a  negative  position),  the  pressurized  section  (starts  at  the 
0-m  position)  is  60  m  long,  and  the  guide  section  is  40  m  long.  Thus,  the  free 
piston  is  initially  at  the  60-m  position.  The  pressure  spikes  on  the  projectile  front 
are  due  to  the  reflection  of  the  traveling  shock  wave  that  borders  between  the 
computational  zones  1  and  2.  In  order  for  the  projectile  to  experience 
deceleration  lower  than  the  allowed  1,600  g,  the  frontal  pressure  on  the  projectile 
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should  not  exceed  38  MPa.  The  free  piston  is  accelerated  rapidly  to  supersonic 
velocity,  but  the  velocity  eventually  decreases  because  of  the  momentum  transfer 
to  the  water  swept  by  the  free  piston.  As  evident  from  Figure  10,  if  the  transition 
tube/rupture  disc  process  is  not  considered  and  the  projectile  is  assumed  to 
impact  directly  on  the  pressurized  gas,  the  eventual  effect  on  the  projectile 
velocity  is  very  small  but  the  projectile  would  experience  significant  higher 
frontal  pressure  peak.  In  all  following  simulations  the  rupture  disc  process  is 
taken  into  account. 


4.3  Effects  of  the  Water  Friction  (Figure  11) 


Figure  11.  Dynamics  in  the  SRS — effects  of  the  water  column  friction. 

Figure  11  demonstrates  the  effect  of  the  water  friction  on  the  projectile  velocity.  If 
water  friction  is  taken  into  account,  the  simulation  indicates  that  the  projectile 
will  exit  the  tube  with  velocity  that  is  about  20  m/s  smaller  than  if  water  friction 
is  neglected. 

4.4  Flow  Dynamics  for  Various  Length  and  Pressure  Parameters 

In  each  of  the  following  simulations,  depicted  in  Figures  12-15,  certain  tube 
lengths  are  selected,  as  well  as  two  initial  tube  pressures  (200  psi  and  300  psi), 
and  the  water  mass  necessary  for  a  projectile  exit  velocity  of  10  m/s  is  foimd 
using  the  Regula  Falsi  [17]  iterative  technique.  In  Figure  12's  caption,  the 
notation  100-60  stands  for  the  total  length  (100  m)  of  the  pressurized  section  and 
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the  guide  section  of  which  the  pressurized  section  (charged  to  either  200  or 
300  psi)  is  60  m.  This  is  the  same  for  the  notations  70-50,  80-50,  80-60,  in  the 
captions  of  Figures  13, 14,  and  15,  respectively.  The  purpose  of  these  simulations 
is  to  find  optimal  tube  lengths  and  pressures.  It  is  desirable  to  have  short  tubes 
that  are  lightly  pressurized.  However,  it  is  not  desirable  to  have  multiple  high- 
magnitude  pressure  spikes  on  the  front  of  the  projectile.  Thus,  a  compromise  has 
to  be  found  with  respect  to  tube  size  and  initial  pressures  and  the  pressure 
loading  on  the  projectile.  Table  1  summarizes  some  of  die  results.  It  is  apparent 
from  Figures  12-15  and  Table  1  that  the  80-60  tube  charged  to  300  psi  is  a  good 
choice.  In  diis  case,  the  maximum  deceleration  on  the  projectile  is  less  than 
1,000  g.  Figure  16  depicts  snap  shots  of  the  pressure  and  velocity  distributions 
in  the  (four)  computation  zones  of  the  80-60/ 300-psi  tube. 

For  this  tube,  the  projectile  wdl  attain  a  relatively  constant  velocity  of  10  m/s 
approximately  12  m  before  the  exit.  These  values  of  distance  and  velocity  are 
conservative  in  the  sense  that  they  allow  for  operational  tolerances  obviating  a 
disastrous  rebound  of  the  projectile  into  the  tube.  Note  that  the  maximum 
pressure  loading  on  the  free  piston  is  much  higher  than  on  the  projectile,  but  this 
is  of  little  concern  as  the  free  piston  is  solidly  built  and  it  is  expendable. 


5.  Summary  and  Conclusions 


A  novel  soft  recovery  system  based  on  ballistic  compression  concept  and 
incorporating  a  free  piston  and  water  is  physically  sound.  Such  a  system  can 
recover  controllably,  within  less  than  120  m,  a  155-mm  (46  kg)  projectile  fired  at 
840  m/s.  Because  of  the  assumptions  involved  in  tiie  formulation,  the  numerical 
simulations  are  conservative  with  respect  to  deceleration  limits  and  projectile 
travel  requirement,  and  thus  actual  performance  is  expected  to  be  better.  The 
simulations  indicate  that  the  projectile  can  be  stopped  using  a  pressurized  tube 
length  of  60  m,  and  a  guide  section  of  40  m.  When  using  a  5-kg  free  piston,  and 
26.84  kg  of  water  distributed  in  the  guide  section,  the  60-m  tube  has  to  be 
pressurized  initially  to  300  psi  (2.07  MPa).  In  this  case,  the  maximum 
deceleration  on  the  projectile  is  less  than  1,000  g. 
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PRESSURE,  MPa 


VELOCITY,  mis 


PRESSURE,  MPa  PRESSURE,  MPa 


Table  1.  Selected  ISS-mm  SRS  parameters. 


Tube  Totah  Length  (m)/ 
Pressurized-Section 
Length  (m)/ 
Pressurization 
(psi) 

Peak  Pressure 
on  the  Projectile 
MPa 
(psi) 

Peak  Pressure 
on  the  Free 
Piston 

MPa 

(psi) 

Water  Mass^  to  Slow  the 
46-kg  Projectile  From 

840  m/s  to  10  m/s 
(kg) 

100/60/300 

23.16  (3,358) 

70.61  (10,238) 

22.72 

100/60/200 

25.93  (3,760) 

56.09  (8,133) 

36.98 

80/60/300 

22.62  (3,800) 

70.85  (10,273) 

26.84 

80/60/200 

29.00  (4,205) 

56.23  (8,153) 

43.72 

80/50/300 

27.75  (4,024) 

76.96  (11,159) 

30.84 

80/50/200 

34.03  (4,935) 

59.89  (8,684) 

45.00 

70/50/300 

29.52  (4,280) 

76.96  (11,159) 

37.13 

70/50/200 

36.74  (5,327) 

59.87  (8,681) 

46.94 

a  From  the  diaphragm  to  exit  (the  transition  length  from  the  vent  holes  to  the  diaphragm)  is 
4.5  m. 

^The  free-piston  mass  is  5  kg. 


shockwave 


Figtire  16.  Snap  shots  of  pressure  and  velocity  distribution  in  the  80-60/300-psi  tube. 
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